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Abstract 

We review three different types of viscoelastic surface instabilities: The Rayleigh - Plateau, the Saffman - Taylor and 
the Faraday instability. These instabilities are classical examples of hydrodynamic surface instabilities. The addition 
of a small amount of polymers to pure water can alter its flow behavior drastically and the type of instability may 
change not only quantitatively but also qualitatively. We will show that some of the observed new phenomena can be 
explained by the use of simple rheological models that contain most of the underlying physical mechanisms leading 
to the instability. A quantitative description however is often only possible close to the onset of the instability or 
for weak deviations from Newtonian behavior. A complete theoretical description is still lacking when the system is 
driven far from equilibrium or for fluids with strong non-Newtonian behavior. 
Resume 

Nous discutons dans ce papier trois types d'instabilites interfaciales dans des liquides visco-elastiques : 1'instabilite 
de Rayleigh - Plateau, l'instabilite de Saffman - Taylor et l'instabilite de Faraday. Ce sont toutes les trois des ex- 
emples typiques d'instabilites hydrodynamiques de surface. L' addition d'une faible quantite de polymeres a de l'eau 
pure suffit a lui conferer un comportement fortement non-newtonien. Les instabilites dans ces solutions sont modifiees 
pas seulement de facon quantitative mais aussi dans leur nature. Nous montrons que des modeles rheologiques sim- 
ples peuvent expliquer l'origine des modifications observees sur les instabilites. Une description quantitative n'est en 
general possible uniquement pres du seuil d'instabilite et uniquement pour des deviations faibles d'un comportement 
newtonien. Une description complete loin de l'equilibre et pour des proprietes non-newtoniennes fortes manque a ce 
jour. 



1. Introduction 

In this paper, we review three typical hydrodynamic surface (or more general interfacial) instabilities: the Rayleigh 

- Plateau instability, the Saffman - Taylor instability and the Faraday instability. These instabilities were chosen 
as they are each representative of a more general class of pattern-forming systems and we have investigated each. 
Here we discuss the modifications of these instabilities that occur when using viscoelastic fluids instead of simple 
Newtonian liquids. 

One of the most prominent examples of a surface instability is a water droplet that detaches from a faucet (for 
a review see HI 0]). The primary stages of the detachment of the droplet can be described by the linear Rayleigh 

- Plateau instability: a cylindrical thread is destabilized by surface tension and splits up into individual droplets, 
thus minimizing the energy of the system. In the non-linear regime, the final detachment or pinch-off of the droplet 
represents a so-called finite time singularity and must be described in a very different way, e.g., by use of so-called self 
similarity solutions. The generation of droplets also occurs in many industrial applications where, instead of simple 
liquids like water or oils, more complex liquids such as polymer melts, polymer solutions, surfactants, suspensions, 
etc. are processed. The hydrodynamics of these liquids differ from those of Newtonian liquids; as a result, the 
character of their surface instabilities can be altered dramatically. When a tiny amount of a flexible polymer is added 
to a solvent, pinch-off is strongly delayed and, in the nonlinear regime, a viscoelastic filament is formed between the 
nozzle and the droplet UlLl2||. Interestingly, this contrasts with theoretical predictions obtained 



from a linear analysis, in which an increase in growth rates and an acceleration of the break-off was predicted 0]. 



Preprint submitted to Compte rendues 



September 28, 2009 



Experiments using yield stress fluids II 1 3L 11411 or granular suspensions 11511 also reveal a strong alteration of the 
detachment process. 

A classical pattern-forming system containing a free surface is the Faraday experiment fnill . When a layer of 
liquid is shaken vertically by a sinusoidal force, the flat interface becomes unstable as soon as a critical amplitude 
is exceeded. Standing waves of high symmetry are formed and patt erns of lines, squares, hexagons, quasiperiodic 
patterns, superlattices and oscillons have been reported 1 17. "l8l Il9l I20I, 21]. Energy is continuously fed into the 
system by the driving force and is dissipated by the viscous flow. This leads to a nonequilibrium situation in which 
the pattern selection process is no longer described by thermodynamic laws but only by the nonlinear dynamics of 
the system. For higher driving strengths, the pattern dynamics may become chaotic or weakly turbulent until droplet 
ejection occurs. Againjthe introduction of a complex liquid can fundamentally alter the pattern selection process of 
the system (HlH&ElllllSlIi. 

Another classical example involving, in this case, a moving interface is the so-called Saffmann-Taylor instability 
1 2911 . When a viscous fluid pushes a less viscous fluid in a narrow channel or Hele - Shaw cell, the interface between 
the two fluids becomes unstable and the formation of air fingers occurs. This instability has received much attention 
as an archetype of pattern-forming systems; it belongs to the class of free-boundary problems, in which the growth of 
the structures takes place in a Laplacian field [30]. A linear stability analysis yields a length scale that in combination 
with the geometry of the system also determines the non-linear growth of the viscous fingers observed. The Saffmann 

- Taylor instability has been intensively studied for Newtonian fluids l3~ll [32L l33ll and lately also for non-Newtonian 
fluids J34 , 35 ] such as polymer or surfactant solution s ll36l I37I l38l l39ll foams, pastes, clays and gels 

and granular suspensions or dry granular media I43I l44l l45l l46ll . The Saffmann- Taylor instability is important for 
a number of applications, including oil recovery, coatings, and adhesive debonding fU These applications often 
involve non-Newtonian fluids, e.g., polymer solutions, which have been shown to fundamentally alter the selection 
process of the viscous fingers (see figure[8]l. 

Further examples of hydrodynamic surface instabilities include the Kelvin - Helmholtz instability (the interfacial 
instability of two liquids that are sheared against each other), the Rayleigh-Taylor instability (that of the interface 
between two liquids of different densities with the heavier liquid placed on top of the lighter) , the Benard Marangoni 
(an instability that occurs when a layer of liquid that is heated from below becomes unstable due to differences in the 
surface stresses caused by thermal gradients), the case in which a thin film flowing down an inclined plane becomes 
unstable against wavy distortion I4H1 and the disintegration of a liquid jet into droplets J2]. Some of these systems 
have been investigated for the case of complex liquids. A special case is the so-called shark-skin instability J49L l50h . 
which refers to the appearance of a wavy distortion on a polymeric fiber that is forced through a hole or slit (a "die") 
above a critical speed. This instability is attributed either to an instability at the solid-liquid interface (stick-slip) or to 
a bulk instability 15 111 . 

For analytical and numerical investigations, treatment of the free interface of these systems is still a challenge. 
Despite the fact that these instabilities have been known for a long time, theoretical understanding of the underlying 
mechanisms was only obtained much later. The Faraday Experiment was the first hydrodynamic pattern-forming 
system that has been described scientifically JT^; however, a successful numerical analysis based on the full Navier 

- Stokes equation was performed only recently 15211 . The pinch-off of a water droplet represents a very general 
problem; however, it was solved analytically only in the 1990is 1 5311 . The Saffman - Taylor instability was described 
by Saffman and Taylor in 1958 12911 . but the mechanism of finger selection in this instability has remained a puzzle 
for several decades. The problem was solved numerically in the '80s 15411 but was only much later solved analytically 
lf55 . 5(1 57]. Analytical treatment of the case of viscoelastic liquids remains challenging even today. 

In this paper, we discuss experimental modification of the Rayleigh - Plateau (followed by the pinch-off of the 
droplet), the Saffman - Taylor and the Faraday instabilities when using one type of non-Newtonian fluids: polymer 
solutions. The addition of small amounts of either flexible or rigid polymers to pure water or water- glycerol mixtures 
can dramatically change the properties of the liquid. We show that the non-Newtonian properties of these solutions 
affect the three surface instabilities in very different ways, and further show that some of the observed new phenomena 
can be explained by the use of simple rheological models describing the hydrodynamics of the complex liquids. 
However, a quantitative description is often only possible close to the onset of the instability or for weak deviations 
from Newtonian behavior. A complete theoretical description is still lacking when the system is driven far from 
equilibrium or for strong non-Newtonian behavior. 
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2. The hydrodynamics of complex liquids 



In this chapter, we discuss the hydrodynamics of polymer solutions and characterize the specific solutions that 
were used in the following studies. From a theoretical point of view, polymers represent very complex systems due to 
the large number of degrees of freedom that exist in their structures. Even a single polymer chain can possess millions 
degrees of freedom due to the many independent molecular bonds between monomers that can rotate more or less 
freely. In equilibrium, entropic forces keep the polymer in a mainly coiled state that allows the highest number of 
conformations. A first approach to representing polymer structure is to model the polymer as a series of connected 
beads. If hydrodynamic interaction between the beads is taken into account, one obtains the so-called Zimm model 
1 58[] . This model predicts that a polymer that is deformed relaxes to its coiled equilibrium configuration with a certain 
relaxation time t p . The relaxation time depends strongly on the molecular weight, the flexibility of the polymer and 
the viscosity of the solvent. It is worth pointing out that the polymer is an entropic spring only, and that in typical 
experiments, energies sufficient to affect the bonds are not reached. In dilute solution, polymers do not interact with 
each other. Higher concentrations can be referred to as semidilute as long as effects such as reptation modes, which are 
observed in polymer melts, are not present. The hydrodynamics of polymer solutions is challenging because it is not 
clear a priori how the coarse graining from the microscopic polymer models to a continuum mechanical description 
can be properly performed. Often, macroscopic models built on a heuristic base that can give some insight into the 
specific flow situation are used. A quantity that is often used to measure the importance of elastic effects under flow 
is the so-called Weissenberg number Wi = T p y, with t p the relaxation time of the solution and y the shear rate of 
the flow. It compares in that way the relaxation time of the polymer to the typical time scale of the flow. Only if the 
Weissenberg number is of 0(1) the flow is strong enough to affect the polymers significantly. 

In our experiments, we deal on the one hand with dilute solutions of the polymers xanthane and PEO (polyethylenox- 
ide) and, on the other hand, with more concentrated solutions of PAA (polyacrylic-co-acrylacid). These polymers are 
important for a wide variety of applications. Xanthane is used to control the rheological properties of products such as 
tooth paste, yoghurt, curd and cream; PEO is used to reduce drag in turbulent flows and PAA is used as a shear thick- 
ener in paper coating applications. These solutions show some of the most frequently encountered non-Newtonian 
properties, making them excellent model systems. We now represent the rheology of these solutions and discuss their 
non-Newtonian behavior in detail. 

Dilute or semi-dilute solutions of the rigid polymer xanthane (M w = 3 * 10 6 g/mol) are shear-thinning. FigureQ] 
shows the shear viscosity obtained for solutions of various concentrations of xanthane. With increasing concentration, 
an increased viscosity at low shear rates and increasingly strong shear thinning when y is larger than a critical shear 
rate are observed. On the microscopic level, one can imagine that the rigid polymers align with the flow at a critical 
shear rate and concentration. The decrease in viscosity can be described by phenomenological models as, e. g., the 
power law model r\ = k\y^ n ~ x \ which is valid over a given range of shear rates. The value of n is a measure of the 
strength of the shear-thinning character of the fluid; for the xanthane solutions, this value ranges from n-0.9 for 50 
ppm (parts per million in weight) to n=0.46 for 1000 ppm. For n-1, the Newtonian limit is recovered. No normal 
stresses perpendicular to the shear plane can be detected in the tested range of shear rates. 

For the flexible polymer PEO, the shear viscosity is nearly constant, whereas significant normal stresses can be 
detected (see figure [2j. With increasing polymer concentration, the Newtonian viscosity increases slightly and the 
normal stresses become more pronounced. The normal stress is a stress contribution normal to the shear plane in 
the direction of the velocity gradient. It is, for example, responsible for the Weissenberg effect: if a rod is used to 
stir a Newtonian liquid in a cup, centrifugal forces depress the surface in a parabolic form. If large normal stresses 
are present in a viscoelastic fluid, the inverse can be observed and the liquid climbs the rod, as can be observed by 
stirring the dough. The existence of normal stresses also indicates a strong resistance to elongation; thus, elongational 
viscosities of fluids that undergo normal stress can be orders of magnitude higher than those of Newtonian fluids. 
For the dilute and semi-dilute polymer solutions presented here, this viscosity is very difficult to obtain by classical 
rheology l59ll . 

On a microscopic level, one can imagine coils of the flexible polymer that are stretched by the shearing flow as soon 
as the shear rate overcomes a critical shear rate, given by the Weissenberg number Wi. The flow properties of solutions 
of flexible polymers can be modeled by a microscopic bead and dumbbell model in which a polymer is represented 
by a spring linking two beads. The dumbbells are convected by the flow. Elastic stretching of the dumbbell due to 
viscous drag and Brownian forces leads to additional stresses in the flow. By averaging and integrating, a continuum 
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Figure 1: Left) Viscosity rj as a function of the shear rate y for different concentrations of xanthane, 50 ppm (•), 100 ppm (▼), 500 ppm (♦) and 
1000 ppm (A). Right) Normal stress N\ as a function of the shear rate y for a xanthane solution of lOOOppm (A). 
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Figure 2: (from 17511 ) Viscosity rj (left) and normal stress N\ (right) as a function of the shear rate y for solutions of different concentrations of PEO 
125 ppm (•), 250 ppm (▼), 500 ppm (♦) and 1000 ppm (A). The solid lines represent a fit of our data to the prediction of the Oldroyd-B model. 
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Figure 3: Complex viscosity rf = rf + irf' for a 2000 ppm PAA (M w = 4 * 10 6 g/mol) in 40/60H'r% water glycerol solution. The intersection of 
the dissipative j]' and the elastic part rf' indicate the inverse relaxation time 1 /t„ of the liquid. 
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Figure 4: Left image: A droplet of water falling from a nozzle. Right image: The same experiment after the addition of lOOppm PEO (M w = 
4 * 10 6 g/mol). 



mechanical model can be obtained. The simplest of these models is the Oldroyd-B model, in which Hookean springs 
with infinite extensibility are assumed. A first outcome of the Oldroyd-B model is indeed the existence of a normal 
stress. One finds a constant viscosity and a quadratic increase of the normal stress with shear rate: 77 = const and 
Ni = ^xj 2 - In this model, the viscosity and Wi, the first normal stress coefficient, are independent from the shear rate 
and are functions of the polymer relaxation time t p . On figure [2] one sees that this model fits the experimental data 
reasonably well. An improved version, the so- called FENE-P model, considers a finite extensibility of the springs 
and allows description of larger deviations from Newtonian behavior. When used in this manner, the FENE-P model 
is also valid for higher shear rates or higher elongation rates. 

For concentrated solutions of PAA, one can measure a linear elastic response for low deformation. The macro- 
scopic Maxwell model assumes that the stresses in the liquid are given by the sum of a viscous part that is proportional 
to the shear rate and an elastic part that is proportional to the shear or deformation. The first consequence is that the 
viscosity becomes a complex number, 77* = 77' + 777", which has both a dissipative and an elastic part and one can 
again attribute a characteristic relaxation time, t p ,to the system. The complex viscosity of the solutions as a function 
of frequency is shown in the figure [3] One can see that even if the complex viscosity cannot be directly fitted with a 
simple Maxwell model, a relaxation time of t p - Is can be extracted from the frequency dependency of the complex 
viscosity at the intersection of the n' and 77" curves. Note that the Oldroyd-B model is also referred to as a convected 
Maxwell model, since for a nonlinear, i.e., finite, deformation, the convection of stress must be taken into account. 



3. Rayleigh - Plateau instability and droplet pinch-off 

3.1. Newtonian liquids 

A liquid column with a free surface always disintegrates into smaller droplets because surface tension leads to 
minimization of the surface-to-volume ratio. Rayleigh showed that the size R of the droplets is determined by the 
wave length A of the sinusoidal distortion with the fastest growth rate. For the inertia-dominated case, i.e., for low 
viscosity liquids, one finds A & 9r with the radius r of the column. For the pending droplet, a similar scenario holds. 
When a pending droplet is fed quasitatically via a syringe pump, it starts to fall if gravitation overcomes capillary 
forces. However, as soon as the droplet beg ins to fall, surface tension again leads to minimization of its surface 
area and acts as the main pinching force |60[ 61, hJ. Thus, the primary stages of the pinch-off of a droplet can be 
modeled as a Rayleigh instability. The thinning dynamics of the neck radius, shown in figure |U can be fitted with an 
exponential law corresponding to exponential growth of the amplitude of the most unstable wavelength. The growth 
rates fit well with the predictions made by Rayleigh's theory ll 
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Figure 5: Left) The minimum neck diameter h m j„ vs. time for water (o) and lOOppm PEO(M„. = 4 * 10 6 g/mol) in water (□). Dotted line: 2/3 
power law line as a guide for the eye. Dashed line: exponential fit corresponding to the Rayleigh-Plateau instability of a liquid column of water. 
Full line: exponential fit in the regime of strongly increasing elongational viscosity. Right) The elongational viscosity for different concentrations 
of PEO as a function of the Hencky strain. 



During the final stages of the pinch-off process in the nonlinear regime, the situation becomes very different. The 
system no longer reflects the original geometry, e.g., the diameter of the nozzle, but instead depends on material 
parameters (density, surface tension and viscosity) only. The dynamics of this finite time singularity in which the 
minimum neck diameter reaches zero in finite time can be described by self-similarity solutions. For the low viscosity 
and thus inertia-dominated regime, which applies in the case of water, it follows that the minimum neck radius 
approaches zero such that r(t) = 0.7 'it - t c ) 2 ^ (figure|Hl i53ll . 

3.2. Polymer solutions 

The effect of elasticity on the (linear) dynamics of the Rayleigh-Plateau instability can be best studied theoretically 
by use of the (linear) Maxwell model. Chang et al. |62l| performed a linear stability analysis and predicted slight 
variations of the critical wave numbers and of the growth rates due to linear elasticity. 

In the experiments presented here (fig§Jl, solutions of the flexible polymer PEO at low concentrations (10-2000 
ppm) in a low viscosity solvent were used. In these solutions, elastic contributions at small deformations are weak and, 
as a consequence, the complex viscosity 77* is difficult to determine by small amplitude oscillatory shear rheometry. 
This means that, for small deformation, the polymers do not affect the flow. Indeed, the experimental results (fi^5]l 
show that the dynamics of the primary stages of the detachment process are not altered. Only when the flow (and 
thus the elongation) is strong enough do the polymers become stretched, interrupting the finite time singularity of the 
pinch-off process. Instead, one observes the formation of a filament and an abrupt transition to a new exponential 
regime with a much larger time scale. This inhibition of the finite time singularity was first observed by Amarouchene 
etal. H. 

This is a surprising observation because it means that the elastic stresses in the liquid that balance the stresses from 
the surface tension are much higher than for the case of the pure solvent (fig. [5}. Still, the shear viscosity of the sample 
is close to the solvent viscosity, even at high shear rates. The solution to this apparent contradiction is the following: 
any type of flow can be divided into a rotational and an elongational part. The elongational part stretches the polymers 
and induces stresses. The rotational part, e.g., in shear flow, causes the polymers to tumble and stresses are averaged 
out to a large degree. In the filament, the flow is purely elongational and stretches the polymers most efficiently. This 
leads to the so-called elongational viscosity r\ e . For Newtonian liquids, the elongational viscosity T] e is directly given 
by geometrical considerations and it follows r\ e - 3Tj s i war . The factor 3 is called the Trouton ratio. For polymers, the 
situation is less clear and measurements of the elongational viscosity are non-trivial. The experiment now gives us 
the opportunity to estimate r\ e . In fact, a similar method is used in a commercialized version of the experiment, the 
so-called CaBER (Capillary Break up Extensional Rheology, Thermo-Fisher Scientific, Karlsruhe, Germany). 

Instead of a falling droplet, CaBER uses a capillary bridge formed between two plates that are separated abruptly 
until the capillary bridge starts to contract. The force balance of the surface stresses and viscous stresses yields an 
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apparent elongational viscosity T] e = ^r, with cr the surface tension 16311 . Indeed, if the filament is cylindrical, the 

elongational rate can be obtained via e = 4| and it is constant for the exponential decrease of the radius found in this 
regime. The minimal neck radius r(t) shrinks exponentially with time and the elongational viscosity r\ e thus increases 
exponentially with time. This is a type of long-lasting start-up situation in which the increase in stresses holds until 
the polymers eventually become fully stretched. Fig. [5]shows the elongational viscosity as a function of the Hencky 
strain H. The Hencky strain, H, is a measure of the elongation and is given by H — J edt. The elongational viscosity 
of the polymer solutions varies from 1 to 100 Pas and is thus five orders of magnitude higher than the shear viscosity 
(see figure[2]l. 

The growing elastic stresses in the filament stabilize the flow, since any distortions would lead to further stresses. 
This makes the filament very robust and we might expect a Rayleigh-Plateau-like instability to be observed only at 
the very end of the thinning process when the polymers are fully stretched and elastic stresses cannot increase further. 
However, in most of the experiments, a more or less irregular instability scenario is observed at the very end of the 
thinning process and singular "beads" grow on the filament. In a numerical study, Chang et al |)62j] predicted that 
the filament should start to disintegrate from both ends where it is connected to the falling droplet and the reservoir 
in the nozzle. At these points, the curvature and surface stresses are maximal. These authors predict an iterative 
disi nteg ration process, similar to the one observed in more viscous Newtonian liquids. Experimentally, Oliveira et 



al. 164 16511 found an iterative process between generations of larger and smaller beads in a CaBER setup, but their 



range of observations was limited. Experimental data obtained by Sattler et al. (66] indeed revealed an instability 
process that was triggered from the ends; instead of an iterative process, we found a Rayleigh-Plateau like instability 
scenario on the viscoelastic filament. In contrast to previous studies, in which the plates where pulled apart abruptly, 
in these experiments the separation was performed very gently in order to prevent any additional distortions. Though 
in most of the experiments singular droplets grew on the filament, it was possible to observe the exponential growth of 
a coherent sinusoidal pattern, (fig[6]right). The inverse growth rate of the pattern was found to be l/co = 9.3 + 0.1ms. 
Linear stability of a viscous fluid thread JlJ] predicts a> = y/(6Ror] e ff), resulting in an estimated extensional viscosity 
of 7] e ff = 9Pas + 2, more than one order of magnitude smaller than the extensional viscosity tje (12/j.m) = lOOPas 
estimated above, but nevertheless four orders of magnitude higher than the shear viscosity. The exponential regime 
was followed by the final shrinking of the remaining thread between the beads, which followed a linear law as expected 
for highly viscous liquids lfl H6lll . 

We must still answer the question of how the filament finally breaks. For polymer concentrations greater than 1000 
ppm, it was observed that the liquid filament connecting two beads did not break; instead, a pattern that resembles a 
solid fiber, as shown in Fig|6]right, a), was formed. The beads sit alongside the filament; experiments on fluid drops 
on a fiber [67] show that there must be a. finite contact angle between the drops and the filament for such a symmetry 
breaking to occur. It follows that the thin filament must have formed a (solid) phase different from that of the drops 
ll68ll . This hypothesis could be confirmed by the scanning electron microscopy (SEM) images shown in Fig. [6]right 
b), in which remnants of two droplets connected by a persistent thin thread are seen. Increased magnification (panels 
c and d) allowed us to estimate the diameter of the fiber as 75 - 150«ra. Our interpretation of this data is that, 
due to the coupling of stress fluctuations and concentration fluctuations, a flow-driven phase separation takes place. 
Microscopically, one can imagine that while solvent drains from the filament, the polymers become entangled, leading 
to even higher polymer concentration and increased entanglement, i.e., flow-induced phase separation takes place ll69tl . 
Further evidence for this concentration process was found in IT70ll . where birefringence measurements were performed 
to examine molecular conformations in the break-up process. Evaporation can be excluded as a factor in the formation 
of solid fibers, based on experiments in a two-fluid system. 



4. Saffmann - Taylor 

4.1. The classical Saffman - Taylor instability 

The classical Saffman-Taylor instability occurs when, for example, air pushes a viscous fluid in a narrow channel 
of height b and width W, a so called Hele-Shaw cell. In the following, the viscosity of air is neglected and the 
viscosity of the viscous liquid is given by tj. The surface tension between the two fluids is cr and the viscous liquid is 
considered to perfectly wet the channel. 
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Figure 6: Left) (Fig. 2 from [66]) Growth of a sinusoidal instability of the viscoelastic filament of an aqueous PEO solution (1000 ppm) that 
develops into a group of droplets on the thinning filament. The spacing of the pictures is 1/300 s. Right) (Fig. 4 from (6fj\ ) a) The final state 
of the filament. Beads are formed off-center relative to the thread, b) Scanning electron microscopy image of two beads connected by a thread 
(intermediate resolution). The structure was caught and dried upon the substrate, c) Another example of the structure; the red box indicates a 
closeup at high magnification shown in d). The diameter of the fiber can be as small as lOnm. 
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Figure 7: Left) Growth rate r as a function of the wave number k. The most unstable wavelength l c is given by 2nlk c . Inset: Snapshots of the 
destabilization of a planar front between air and silicon oil (courtesy D. Derks and A. Lindner). Right) relative finger width A as a function of the 
control parameter 1/B for air pushing silicon oils of different viscosities in channels of different geometries. Inset: snapshot of a finger advancing 
into a linear cell (courtesy C. Chevalier). 
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Flow in the confined geometry is then governed by Darcy's law, which gives the mean velocity (averaged over the 
thickness of the channel) of the fluid as a function of an applied pressure gradient : V = -j^-Vp- The incompressibil- 
ity of the fluid reads V ■ V = and one thus deals with growth in a Laplacian pressure field S.p = 0. The pressure jump 
at the interface is given by 5p = cr (| + k), with k being the curvature in the direction of the channel width, once again 
using a two-dimensional approximation. Together with the boundary conditions, this set of equations completely 
determines the problem. 

When the less viscous fluid pushes the more viscous fluid, an initially straight interface becomes unstable. Small 
perturbations lead to an increased pressure gradient and a higher velocity in front of the perturbations and are thus 
amplified. Surface tension, on the other hand, stabilizes the initially straight interface. The competition between 
viscous and capillary forces leads to the emergence of a characteristic length scale that can be calculated using linear 
stability analysis 17111 . The maximum growth rate is found for a wavelength l c = nbj 'sfCa with capillary number 
Ca - rjU/cr. The small fingers grow and begin to compete with the more advanced fingers, screening the less 
advanced fingers. In a linear channel of width W, one finally observes a single finger propagating through the cell, 
the result of a non-linear growth process. An example of initial finger growth and finger competition, together with 
the growth rate obtained from the linear stability analysis, can be seen in figure [7^. The relative width of the single 
finger A, defined as the ratio between the finger width w and the cell width W, is given by the control parameter 
l/B = l2Ca(W/b) 2 ; that is, the ratio between the two length scales of the system W and l c . Representation of the 
results obtained using different fluids (and thus different surface tensions or different viscosities) and different cell 
geometries as a function of the control parameter l/B shows that the results fall on a universal curve (see figure |7J)). 
With increasing velocity U of the finger tip, viscous forces become increasingly important compared to capillary 
forces and the relative finger width decreases. At high velocity, the finger width does not, however, tend to zero but 
stabilizes near a plateau value at A = 0.5. 



This instability was described by Saffman and Taylor 112911 in 1958; however, finger selection remained a puzzle 



for several decades. Neglecting surface tension, Saffman and Taylor found a family of analytical solutions of the 
shape of the interface given by x = In (l + cos ^l)] that agrees well with experimental observations. This 

treatment does not, however, explain the selection of a given finger width. For this, one must take the surface tension 



into account. This was done numerically by McLean and Saffman in 1981 115411 . The selection process was solved 
analytically only much later 115 5l l56l 15711 and was attributed to the fact that the surface tension represents a singular 
perturbation leading to a solvability condition at the finger tip. It is this condition that selects the finger from the 
family of solutions found by Saffman and Taylor. 

So far, we have treated an ideal two-dimensional problem. In reality, however, there is a thin wetting film that 
remains between the advancing finger and the glass plates. The thickness of this film, according to the Bretherton 
law, is proportional to Ca 2 ^ t/R = 0.643(3Ca) 2 ^ 3 17211 . As a consequence, the pressure jump at the interface is 
continuously modified. This three-dimensional effect leads to the slight modification of the finger width observed 
between different experimental geometries, as observed in figure [7] 

4.2. Finger narrowing in shear thinning fluids 

When performing experiments in dilute solutions of xanthane, one observes a strong modification of the selection 
process. As can be seen in the snapshot in Figure^, at high velocity, fingers are found to be significantly narrower 
than the classical limit of A — 0.5. This observation can be qualitatively linked to the behavior of the shear thinning 
fluid pushed by the finger in the Hele-Shaw cell. Numerical simulations Jz3 , 74 ] show that the viscosity is not uniform 



throughout the cell (see figure[8j)); regions of high fluid velocity and thus high shear rate have a low viscosity. This is 
essentially the case in front of the finger tip and the system becomes anisotropic, leading to finger narrowing. 



For weak shear thinning (1 > n > 0.65), we have shown 175L 13611 that simply replacing in the control parameter 
l/B the constant viscosity rj by a shear- dependent viscosity rj{y) allows rescaling the data onto the universal curve for 
Newtonian fluids. The shear rate y is here the average shear rate in the cell. For stronger shear-thinning (n < 0.65), 
this rescaling fails and deviations from the classical result toward smaller fingers are observed. 

Narrower fingers have also been observed by Rabaud et al. lT76Tl . Using a Hele - Shaw cell with engraved glass 
plates, they found viscous fingers with A significantly smaller than 0.5 for Newtonian fluids. The observation of such 
"anomalous" fingers is explained by the fact that the engravings represent a local perturbation at the finger tip. This 
disturbance removes the classical selection of the discrete set of solutions. The continuum of solutions given by the 
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Figure 8: Left) Relative finger width A as a function of the velocity v for solutions of xanthane (o) at 1000 ppm and PEO (•) at 50 ppm. Inset: 
snapshots of fingers in the two solutions at high velocity (4 cm/s). Right) Viscosity in front of the advancing finger in a shear thinning fluid found 
from numerical simulations (73ll . The shear thinning effect increases from top to bottom. Red corresponds to high viscosities and yellow to low 
viscosities. 

analytical result of Saffman and Taylor without surface tension then becomes accessible: A can take values smaller 
than 0.5 at high velocity. Rabaud et al. showed that for a given value of the capillary number Ca, it is not the relative 
finger width that is selected but that the dimensionless radius of curvature at the tip p/b. p can be linked to the finger 
width A via the relation p = , which follows from the finger shape predicted by Saffman and Taylor. 

A similar mechanism has been found to be responsible for the selection of the viscous fingers in a shear thinning 
fluid. Here when the shear thinning character of the fluid is strong enough, anisotropy plays the role of the perturbation 
at the finger tip. Figure [9^ shows experimental finger profiles for a given capillary number and three different cell 
widths. One clearly observes that the radius at the finger tip is identical for the three experiments, leading to lower 
finger width compared to the Newtonian case, in which the relative finger width A is selected. The relation between 
p/b and A is found to depend on the shear thinning character of the fluids and thus the shear thinning exponent n. 
Knowledge of the relationship between p/b and A solves the selection problem, as we can now predict the finger 
width A from the rheological data. The presence of shear thinning thus leads to a completely different selection 
mechanism that is closer to what is observed, for example, in dendritic growth lf3Qh and which requires anisotropy 
in the system. Corvera Poire et. al. directly solved the problem for a power law fluid; their results are in good 
agreement with the experimental observations. Note that at high velocities one observes a saturation of p, leading to 
an increase of the finger widths. This might be attributed to inertial effects, which have been observed to increase 
finger width [43] and which may begin to play a role in this low viscosity fluid. 



4.3. Finger widening due to normal stresses 

In experiments using solutions of the flexible polymer PEO, completely different behavior is found. In contrast to 
the observations in shear thinning fluids, where finger narrowing occurs, one now observes finger widening compared 
to the Newtonian case (see figure^). 

The presence of normal stresses in the thin wetting layer might be responsible for the finger widening; one can 
attempt to account for this effect by adding a supplementary pressure to the system. In classical theory, the pressure 
jump at the interface between two liquids is given by the radius of curvature. Tabeling et al. 17211 have shown that 
one can incorporate the effect of a finite thickness of the wetting film by correcting the surface tension in the control 
parameter. Following the same argument, in the control parameter one can add the supplementary pressure caused 
by the normal stresses to the surface tension term cr* = o~ + \/2N\(y)b. For moderate normal stresses, this allows 
rescaling of the data onto one universal curve and once again solves the selection problem. 
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Figure 9: Left) Experimental ringer shapes in a solution of xanthane of 2000 ppm for three different channel widths W=2, 4, and 8 cm. Right) 
Sketch of the thin wetting layer observed between the finger and the glass plates 



5. Faraday instability 

5.1. History 

The Faraday experiment was first reported in 1831 l[l6ll . In an appendix of a paper on Chladni figures, Faraday 
reports on the crispated state of a layer of liquid that is shaken vertically. He refers to works by "Oersted, Wheatstone 
and Weber and probably others (sic)" who had earlier mentioned the phenomenon; however, according to Faraday, 
it was he who gave the first conclusive description. The experimental setup consisted of a box that was mounted 
on a rod. The rod was set into vibration by a bow; its oscillation frequency was presumably on the order of a few 
Hz. When the oscillation amplitude exceeded a certain critical value, standing capillary surface waves were observed. 
Apparently, Faraday was impressed by the richness of the patterns and noted: "obtained in this way the appearances 
were very beautiful, and the facilities very great". It is a remarkable achievement that even at that early date Faraday 
found that the surface waves oscillate with half of the driving frequency. This is the so-called subharmonic response. 
Faraday also pointed out differences in the wavy surface patterns when he compared simple oils with, e.g., the white 
of an egg. He stated: "The difference between oil and white of egg is remarkable; . . . the crispated state may be a 
useful and even important indication of the internal constitution of different fluids." 



In 1868, Matthiesen [78] reported on systematic measurements; he stated incorrectly that the surface response 



should be synchronous to the drive. In 1883, Lord Rayleigh [79] proposed a theoretical treatment in terms of a 
parametric pendulum, the Mathieu oscillator. In 1954, Benjamin and Ursell lioll solved the linear problem for ideal 
liquids (without viscosity) with an infinite set of Mathieu oscillators that oscillate with integral and (« + l)/2 multiples 
of the driving frequency. The integral multiples correspond to a possible harmonic (synchronous) response and the 



(n + l)/2 to a subharmonic response (figurefTOIi. In 1994, Kumar and Tuckermann 118 111 presented a numerical analysis 
of the linear problem in the case of viscous liquids with finite depth of the layer. Together with an analytical treatment, 
this analysis was used to find parameters to experimentally reproduce a harmonic response using a very thin layer of 
liquid 18211 . Recently, a large variety of patterns with up to 10-fold rotational symmetry (quasiperiodic), superlattices 
or localized patterns have been reported M20l Il9l 18311 . A variety of these patterns could be obtained using a simple 
liquid driven with a single frequency. Near onset, in linear order, a single wave number first became unstable. The 
resulting wave vector(s) could have any orientation, but nonlinear interaction with the higher harmonics led to a given 
pattern selection process II 17[l82l 18411 . Obviously, the pattern dynamic might become even richer if two discrete wave 
numbers become unstable simultaneously. 

It is worth mentioning that numerical simulation of the full hydrodynamical problem of the Faraday experiment 



with a single driving frequency for simple liquids became available only very recently [52]. Perinet et al. solved the 
complete nonlinear Navier - Stokes equations by a finite-difference projection method coupled to a Front Tracking 
technique for the calculation of the surface tension forces and advection of the interface. They compared the complete 
spatial and temporal Fourier spectrum of the surface state and found good agreement with experimental data from 
Kityk et al. l85ll . A quantitative theoretical description of the nonlinear wave state of complex liquids remains an even 
more challenging problem that is still unsolved. 
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Figure 10: (From (27)]) The linear stability diagram of the Faraday experiment for a Maxwell fluid. Full lines mark the transition to the subharmonic 
response and dashed lines to the harmonic one. a) For a Newtonian fluid, i.e., when the polymer relaxation time t p is zero, the subharmonic response 
becomes unstable first, b) For finite relaxation times, a new harmonic tongue appears c) Only when the inverse polymer relaxation time compares 
to the driving angular frequency CI does the new harmonic instability tongue become unstable first. 



5.2. Complex liquids 

The first experimental data on Faraday waves of polymer solutions were, to our knowledge, presented in 1998 by 
Raynal, Kumar and Fauve ll23ll . Their work concentrates on dilute polymer solutions where the influence of elasticity 
is small. They found a slight shift of the critical acceleration; the critical wave numbers were not affected. In 1999, 
Miiller and Zimmermann presented a linear stability analysis for a Maxwell fluid 112711 . They found that when the 
inverse of the relaxation time of the Maxwell fluid compares to the driving frequency, a harmonic response might 
become unstable first and that, for a certain set of parameters, a bistable situation exists; see fig. \TU\ Again, in 1999, 
Wagner and Miiller presented experimental and theoretical data on a Faraday experiment using a concentrated poly- 
mer solution of 2000 ppm PAA (M w = 4 * 10 6 g/mol) in a water-glycerol mixture 112411 with a relaxation time t p ~ 1 s 
(see section [2j. The rheological data obtained for the complex viscosity was then fed into the linear stability algo- 
rithm of Kumar and Tuckerman 118 ill and a reasonable agreement of the experimentally and numerically determined 
critical accelerations was obtained. Notably, a bicritical situation was found in which the subharmonic and harmonic 
responses became unstable simultaneously. 

In fig. [TT] an overview of the observed patterns is presented. This phase diagram was obtained by setting a fixed 
frequency and varying the driving strengths. When a (new) pattern was visible, the transition point was noted. The 
primary instability was found to be hysteretic (e < 0) and most of the transitions between the patterns were hysteretic 
as well (not shown). For frequencies well above the bicritical frequency /j » 39Hz, a subharmonic pattern of lines 
(region V) was observed corresponding to observations in a Newtonian liquid of similar viscosity. For parameters 
where only a harmonic response existed, a coherent pattern of hexagons was observed (region Ha, fig. [T2"h). The 
hexagonal symmetry is generic for a harmonic response where the temporal symmetry allows the coupling of three 
wave vectors H„ that are equally distributed on the critical circle (figfT2b). For frequencies close to the bicritical point, 
with increasing driving strength subharmonic wave vectors S „ also become unstable. They are slaved by the harmonic 
hexagonal pattern and arrange together to a 2 x 2 superlattice. The nomenclature is taken from crystallography and 
relates to the ratio of 1 : 2 of the harmonic and subharmonic wave vectors, which is close to the ratio of linear unstable 
wave numbers (region Ilia, fig. [12b). 

The pattern-forming process can be understood in the following way. The linear stability analysis reveals that 
the surface state consists of a singular wave number k (only exactly at the bicritical point do two wave modes syn- 
chronously become unstable) but of an infinite series of temporal Fourier components, «Q or (« + 1)/2Q for the 
harmonic or subharmonic responses, respectively. This linear result is fed into the nonlinear equations for the hydro- 
dynamic velocity field v. The solvability condition in quadratic order implies that resonant terms must vanish in the 
higher orders to prevent secular growth. For the temporal components, an arbitrary quadratic nonlinearity results in 



12 



-5? 

03 
i 

03 



0.10 
0.08 
0.06 
,0.04 
,0.02 
0.00 



II 

w -0.02 
-0.04 
-0.06 



. IVa 

"\ 

- \ 

' E 

\ Ilia 
\ 

A— 


IVb 

\ 

□ 


\ 

IVc 


i i i 
♦ — » — 


" I la V 


lllb\ 
lb 


_ - - - * 

□ 

y 

/ 

-□ 

i i i i i 



30 35 40 45 50 55 60 

forcing frequency f=Q/2n (Hz) 



65 



Figure 11: (Fig. 2 from |24|| ) Phase diagram of the observed nonlinear patterns. The symbols mark experimental data points. The abscissa, 
e = ^--^ = o, indicates the linear threshold with a the acceleration amplitude and a c the critical acceleration. The lowest solid line denotes the 
saddle point of the hysteresis. Region I: flat surface; Ha: harmonic hexagons covering the whole surface; Ilia: harmonic-subharmonic hexagonal 
superlattice (see text) extending over the whole surface; lib, Illb: as before, but patterns occur in the form of localized patches surrounded by the 
flat surface; IVa, IVb, IVc: chaotic dynamics of subharmonic lines competing with extended (a) or localized (b) hexagonal superlattices or the flat 
surface (c); V: stationary subharmonic lines extending over the whole surface (Newtonian regime). 




Figure 12: (a) The harmonic hexagonal pattern, (b) the harmonic-subharmonic hexagonal superlattice. (c) and (d) the respective spatial Fourier 
spectra, (e) (Fig. 4 from (24}l) Localized stationary surface patterns of harmonic hexagons in coexistence with a localized nonstationary patch of 
lines. 
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a frequency spectrum of integral multiples of f2, whether or not S or H are considered. Thus, quadratic nonlineari- 
ties are able to resonate with harmonic linear eigenmodes, but not with subharmonic ones. In the same way, spatial 
resonance must be guarantied as well. Now, any triplet of harmonic modes {^HX^Hi^m) with |k# m | = kn and 
k//i + k#2 + k//3 = are in resonance. This generic 3-wave vector coupling is well known (e.g., from non-Boussinesq- 
Rayleigh-Benard convection) and enforces a saddle node bifurcation towards hexagonal patterns. The associated 
solvability condition is referred to as the amplitude equation or Ginzburg-Landau equation. Within the subharmonic 
regime, such a resonant 3-wave vector coupling is prohibited due to the missing temporal resonance. The pattern se- 
lection mechanism is therefore then controlled by the cubic coupling coefficient in the associated amplitude equations 
1184 18611 . and a variety of patterns is allowed. 

The secondary superstructure must result from a nonlinear excitation process, as well. The responsible mechanism 
is again a 3-wave vector interaction. The linear analysis yields the result that the wave number of the subharmonic 
response is approximately |ksi| = ^|k//i|. Invoking basic resonance arguments, one finds directly k,si = ^k^i- 

For increased driving strength, localized patches of hexagons have also been observed. They cannot be explained 
by a triplet of real Landau equations supplemented by diffusive spatial derivatives, since within this familiar amplitude 
equation model, stable isolated islands of hexagons do not exist. For even larger driving strength, the patterns become 
chaotically time-dependent. Patches of subharmonic ally oscillating lines originating in an erratic manner from the 
cell boundary or the flat surface penetrate into the stationary hexagonal superlattice. They then disappear and the 
original structure is recovered. This process repeats itself on time scales of seconds to minutes, leading to a temporary 
coexistence of the stationary hexagonal superlattice with subharmonic lines (region IVb, fig. \12[ . Higher driving 
amplitudes lead to a fully chaotic surface pattern. Still, no satisfying theoretical description of these states exists. 



6. Summary 

In the present paper, we have revisited three classical surface instabilities, the Rayleigh-Plateau instability fol- 
lowed by the detachment of a droplet, the Saffman-Taylor instability and the Faraday experiment. We have described 
the modifications of the selection processes that occur when polymer solutions instead of simple Newtonian fluids 
are used. We have shown that the non-Newtonian properties that come into play when using viscoelastic fluids can 
fundamentally modify the selection processes. The Rayleigh-Plateau instability, which describes the onset of the 
detachment of a droplet, is not modified by the use of dilute solutions of a flexible polymer, as the linear viscoelas- 
tic properties of this solution are weak. The final detachment of the droplet is, however, strongly delayed as high 
stresses develop when the fluid filament becomes more and more stretched. The additional resistance to breakup can 
be attributed to large elongational viscosities present in these solutions. Only in the final stages of the detachment 
process, when the polymers are fully stretched, might another Rayleigh-Plateau instability be observed. Eventually, 
a flow-induced phase separation may lead to a solid nanofiber. Using this type of polymer solution also leads to a 
modification of the Saffman-Taylor or viscous fingering instability. Normal stresses develop in the thin wetting film 
left between the advancing finger and the glass plates of a Hele-Shaw cell, leading to finger widening compared to 
the classical results. When using shear thinning solutions of a rigid polymer, finger narrowing is observed. This can 
be attributed to an anisotropy that develops in the system. In front of the finger tip, where the shear rate is high, the 
viscosity is decreased. For strong shear thinning, the selection process is quantitatively modified as the anisotropy 
removes the solvability condition at the finger tip that is responsible for the selection process in Newtonian fluids. The 
selection process resulting from this anisotropy is found to be closer to that involved in the formation of dendrites. The 
Faraday instability is studied in a concentrated polymer solution that can be approximately described by a Maxwell 
model. It has been shown that when the relaxation time of the complex fluid is of the same order of magnitude as the 
driving frequency, a subharmonic and an harmonic response can become unstable simultaneously. This leads to an 
entirely new class of patterns, the so-called superlattices. 
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